ISAS/EP/98/63 



Exact Matrix Elements in Supersymmetric Theories^ 



G. Mussardo 

International School for Advanced Studies, Trieste 

and 

Laboratoire de Physique Theorique et Hautes Energies, Paris VI 



Abstract 

The lowest representatives of the Form Factors relative to the trace operators of 
N = 1 Super Sinh-Gordon Model are exactly calculated. The novelty of their de- 
termination consists in solving a coupled set of unitarity and crossing equations. 
Analytic continuations of the Form Factors as functions of the coupling constant 
allows the study of interesting models in a uniform way, among these the latest 
model of the Roaming Series and the minimal supersymmetric models as investi- 
gated by Schoutens. A fermionic version of the c-theorem is also proved and the 
corresponding sum-rule derived. 
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1 Introduction 



The solution of a Quantum Field Theory (QFT) is provided by the complete set of correla- 
tion functions of its fields. Perturbation theory often proves to be an inadequate approach 
to this problem therefore more effective and powerful methods need to be developed. In 
this respect, one of the most promising methods is the Form Factor approach which is 
applicable to integrable models [p], 0. This consists in computing exactly all matrix ele- 
ments of the quantum fields and then using them to obtain the spectral representations of 
the correlators. In addition to the rich and interesting mathematical structure presented 
by the Form Factors themselves (which has been investigated in a series of papers, among 
which [2-13], the resulting spectral series usually show a remarkable convergent behaviour 
which allows approximation of the correlators (or quantities related to them) within any 
desired accuracy 0, 0, 0, 0, |T|, |]|. 

Particularly important examples of QFT are those which are invariant under a super- 
symmetry transformation which mixes the elementary bosonic and fermionic excitations. 
The subject of this paper is the investigation of the Form Factors of bidimensional su- 
persymmetric theories, in particular of the Super Sinh-Gordon Theory (SShG) and of 
those models which can be obtained by its analytic continuations. In these models, the 
degeneracy of the spectrum dictated by the supersymmetry implies the existence of multi- 
channel scattering processes and the resulting S"-matrix is necessarily non-diagonal. In 
this case, the complete determination of the matrix elements of the quantum fields for an 
arbitrary number of asymptotic particles becomes a mathematical problem of formidable 
complexity. Although this problem has not been solved in its full generality in this paper, 
however some progress has been nevertheless achieved in determining the lowest Form 
Factors of these non-diagonal scattering theories. In particular, a set of functions has 
been identified which might be used in future studies to achieve the determination of all 
Form Factors of the models. As we will see, the calculation of the lowest matrix elements 
of supersymmetric theories presents some novelties which makes their determination an 
interesting issue in itself. In fact, for the first time one has to solve a coupled set of 
unitarity and crossing equations which originate from the fermionic nature of the super- 
symmetric theories. It should be said that in the past the calculation of the lowest Form 
Factors of supersymmetric theories has been approached in a paper by Ahn |24] but it 
turns out that his results were incorrect, as briefly discussed in Appendix B. 

The paper is organised as follows. The next section recalls the basic results of the 
scattering theories of bidimensional models with a iV = 1 supersymmetry. Section 3 deals 
with the simplest super-symmetric scattering theory which will be finally identified with 
a particular point of the ordinary Sine-Gordon S'-matrix. The Super Sinh-Gordon Theory 
(SShG) and the relevant features of the deformed superconformal models are the subject 
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of Section 4. The S'-matrix of the SShG model and of others which can be obtained as 
analytic continuations thereof are analysed in Section 5. The determination of the lowest 
Form Factors of a particularly important class of operators - those of the trace of the 
super stress-energy tensor - is discussed in Section 6. In Section 7, a thorough check of 
their validity is obtained by means of the c-theorem sum rule as well as by comparing 
them with the limiting form obtained in two significant cases, the lastest representative of 
the Roaming Models and the first example of the Schoutens's models. Section 8 is devoted 
to the fermionic formulation of the c-theorem which can be achieved in supersymmetric 
theories whereas Section 9 contains the summary of the results and conclusions. There 
are two appendices, the first relates to the properties of a function entering the Form 
Factor calculation, the second contains a brief discussion and criticism of previous results 
obtained by Ahn. 



2 Generalities and Notation 

Scattering theories of integrable super-symmetric theories have been discussed in detail 



in Schoutens's paper |20| (see also pl|). However, for the sake of clarity, in this section 
we explicitly state all notations and conventions which will be used in this paper. 

Let us consider a two-dimensional quantum field theory made of a bosonic and fermion 
particle both of mass m. The one particle-state of the bosonic and fermionic particles will 
be denoted by | b{j3)) and | /(/?)) respectively^], where (3 is the rapidity, i.e. the variable 
entering the dispersion relations p = m cosh (3 and pi = m sinh (3 (in the following we 
will also consider the combinations p± = po ± pi = me ±l3 ) . 

Let us assume that such a theory is both integrable (i.e. there exists an infinite 
number of conserved charges) and invariant under a N = 1 supersymmetry. This means 
that among the set of conserved quantities there are two fermionic charges Q and Q which 
satisfy the relations 

Q 2 = P + , Q 2 = P_ , {Q,Q} = , (2.1) 

where P± are the right /left components P± = Pq ± P\ of the momentum operator. The 
action of these charges on the one-particle states can be represented as 

Q | &(/?)> = u v^e^ 2 | /(/?)) ; Q | &(/?)} = 5 ^Ee~^ | /(/?)) ; 
Q | /09)> = c^vW 3 / 2 | b{(3)) ; Q \ /(/?)) = r 1 ^^ 2 I K(3)) , 



1 For the multi-particle states we have | ^(/Ji)^^) ■ • • A n ((3 n )}, where each Ai is either a b or a 
/ particle and the rapidities are ordered in increasing order Pi > fii > ■ ■ • Pn for the in-state and in a 
decreasing order /3i < < ■ • ■ 0n for the out-state. 
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i.e. in terms of two matrices 

Q = 



OJ 



u 
1 



, Q 



5 

5- 1 



(2.3) 



satisying Q 2 = Q 
condition 



1. The anti-commutativity of the operators Q and Q gives the 



(2-4) 



with the actual values of uj and 5 which will be fixed by later considerations. 

The action of Q and Q on a multi-particle states must take into account the fermionic 
nature of these operators and therefore involves brading relations. In the notation of 



ref. ||21|| , we have 

Q | A 1 (p 1 )A 2 (p 2 ) . . . An(p n )) = v^E e&/2 



(2-5) 



k=l 



(Q F A 1 ((3 1 ))(Q F A 2 ({3 2 )) . . . (Q F Ak-i(Pk-i)(QA k (f3 k ))A k+1 ((3 k+1 ) . . . A n (f3 n )) 



and 



Q | A 1 (J3 1 )A 2 (J3 2 ) • • • A n ((3 n )) = V^T, e ' M2 

k=l 

| (Q F A 1 (/3 1 ))(Q F A 2 (p 2 )) . . . (Q F A k _i((3 k -i){QA k (!3 k ))A k+l ((3 k+1 ) 



(2-6) 



MM 



Qi 



(2-7) 



where Q F is the fermion parity operator, which on the basis | b) and | /) is represented 
by the diagonal matrix 

' 1 
-1 

Particularly important is the representation of the two super-charges on the two-particle 
states | b((3i)b((3 2 )), | f(flx)f(fo)), I /(A )&(&)>, I Kfc) f {(3 2 )) . The first two states belong 
to the F = 1 sector whereas the remaining two states to the F — — 1 sector. By choosing 
for convenience /?i = (3/2 and (3 2 = — (3/2, the operator Q will be represented by the 
matrix 



/ 



W) 



o 










LUX 



\ 



-uj 1 x 1 UJ l x 



UJ 1 X ~UJX 1 



\ UJ 1 x 1 



UJX 










/ 



(2-8) 



where x = e^/ 4 . For Q we have analogously 

/ 









dx- 1 


Sx 










d^x- 1 


5- l x- 1 


— 5x 








5- l x 


5x~ l 









(2-9) 
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In the following we will also need the representation matrix of the operator QQ on the 
above two particle states, given by 



/ 



(QQ)(P) 



V 



i sinh | 





-u8 sinh | 












-fcoshf 






-fcoshf 







(2-10) 



Let us consider now the structure of the elastic two-body S-matrix of the theory. 
Since fermions can only created or destroyed in couples, we have the following scattering 
channels 



/0W(&)> = C((3 12 ) I &(&)&(&)> + D(p 12 ) I /(ft)/(/3i)> 

m)w)) = e (m i fmwiV+Fifa) i www) 

b{Pi)f(M = G{Pu)\f(fa)b(fa)) + H{Pu)\b{p2)f{Px)) 



(2.11) 



where (3\ 2 = fi\ — (3 2 - Respecting the ordering of the 2-particle states as above, the 
S-matrix can be then represented by 



( A 


B 










c 


D 
















E 


F 




V o 





G 


H 


) 



(2.12) 



The amplitudes A((3) and D{j3) describe the transmission channels of two bosons and 
two fermions respectively whereas the amplitudes B{(3) and C(f3) are related to the 
annihilation-creation channels of two bosons into two fermions and viceversa. The other 
amplitudes E{(3) and H(f3) describe the reflection channels of the boson-fermion scatter- 
ing whereas the remaining F(/3) and G(/3) describe the pure transmission processes in the 
F — — 1 sector. 

Let us now impose the invariance of the scattering processes under the action of the 
supersymmetric charges. First consider the invariance of the S-matrix with respect to 
the operator QQ: with the choice fti = (3/2 and (3 2 = —(3/2, by taking into account the 
different ordering of the rapidities of the out-states, this implies the following constraint 
on the scattering amplitudes p0| 



QQ((3)S((3) = S((3)(QQ)(-(3) . 
In the F — — 1 sector, this equation leads to the conditions 

E((3) = H((3) , F{p) = G{p) , 



(2.13) 



(2.14) 
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whereas in the F = 1 sector, eq.( |2.13| ) implies 

B((3) = (uj5) 2 C((3) 







2u 2 C{(5) = [A(J3) + D(J3)] sinh | ; 
|s(/3) = [A((3) + D((3)}smh^ . 



(2.15) 



Hence, we can choose B(f3) = C(f3) provided that the constants uj and 5 satisfy the 
condition 



(ujSf 



(2.16) 



Since from the closure of the supersymmetric algebra uj = ±i5, it follows that 5 has to 
satisfy the condition 5 4 = — 1 and therefore we have the following possibilities: 



5 = e^ 



or 



5 



— 7 — 

e 4 



uj = ±e 1 * 



uj = ±e J 4 



(2.17) 
(2.18) 



In the following we will adopt the values uj = e*4 and S = e l ±. With this choice, we have 
then 



B(P) = C(P) , 2B(/3) = -i[A{P) + D{(3)} sinh | . 
In light of all the above equations, the S matrix can be written then as 



S{0) 



(2.19) 



( A 


B 





\ 


B 


D 














E 


F 


V o 





F 


E ) 



(2.20) 



Further conditions on the above functions are provided by the invariance of the S"-matrix 
under each supercharge separately (having already considered the invariance under QQ, 
one can consider only one of them, say Q) in the form 



Q{p)S{p) = S{P)Q{-P) 
This equation is equivalent to the following conditions 



(2.21) 



A{(3) - D{(3) 
A{p)+D{p) 

E{f3) tanh| 



2F((3) 
2 

cosh | 



E((3) 



(2.22) 
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Finally, crossing symmetry implies 



A{in 


~P) 


= A(fi) 


D{m 


~P) 


= D{(5) 


F(l7T 


~P) 


= m 


E(m 


~P) 





whereas the unitarity condition S(/3)S(—/3) = 1 gives 



A(P)A(- 


-(3)+B{(3)B{- 


~P) 


= 1 ; 


A{p)B(- 


-/3) + B(P)D{ 


~P) 


= ; 


B((3)A(- 


-P) + D{(3)B{ 


-P) 


= ; 


B([3)B( 


-(3)+D{(3)D{ 


-P) 


= 1 


E(/3)E(- 


-(3) + F((3)F( 


-P) 


= i ; 


E((3)F(- 


-(3) + F((3)E( 


-P) 


= . 



(2.24) 



In the next section we will discuss several examples of scattering theories which fulfill 
the above set of equations. 



3 The simplest SUSY S-matrix 



The simplest supersymmetric S-matrix can be obtained by noticing that the amplitude 
F(f3) is invariant under crossing. This allows us to make the consistent choice F{(3) = 0. 
In view of the equations ( |2.22j) , the 5* matrix can be written in this CcLSG cLS 



/ -1/coshf 



S(P) 



V 



i tanh f 









i tanh f 
-1/ cosh „ 









-1 




\ 



-1 / 



R(P) 



(3.1) 



where the function R{(5) has to satisfy the equations 



R(P)R(-P) 

R{lTt -P) = 



: tanh f R(P) 



(3.2) 



There is a geometrical method of solving the above coupled set of equations for a mero- 
morphic function with only poles and zeros, none of which are in the physical strip 
< lm(3 < 7r. First of all, the poles and zeros of R(P) are linked to each other by 
the first eq. of ( |3.2|) because if the function R(f3) has a pole (zero) in inr], it should nec- 
essary have a zero (pole) in —mr\. On the other hand, changing the rapidity (3 according 
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to (3 



ITT 



(3 is a reflection with respect to the point in. The function — itanh^ has 



zeros in (3 — ±2kni and poles in (3 
representation is given by 



±(k + l)7ri (A; = 0, 1, . . .) and its infinite product 
[2Jbr-i/3][(2(Jfe + l)7r + i/3] 



L [(2k + l)n + t(3}[(2k + 1)tt - if3] 



(3.3) 



Since the second equation can be written as 

R(in - (3)R(-(3) 



-i tanh 



(3.4) 



to solve this equation we therefore need to find a function whose sovraposition of poles 
and zeros coming from the combination of the two reflections (3 — > in — (3 and (3 — > —f3 
match those of (|3.3|) . It is simple to see that such a function should have zeros and poles 
of increasing multiplicities and its infinite product representation is given by 



R((3) 



n 

fc=0 



' [(2k + 1)tt + 0\ [(2k + 3)n + iJ3] [2(k + 1)tt - i(3} 2 
[(2k + 1)tt - i(3] [(2k + 3)vr - i/3] [2(k + l)n + i(3] 2 

The above function can also be written as 



k+l 



(3.5) 



m = n 



^ T(k + l + i£)r(k + l + i£) p (k + i-i£) 



or as 



R((3) = exp 



1 r°° dt sin 

2 To 



l3t 



t cosh 



2 t 



2 J 



(3.6) 



(3.7) 



It is now interesting to note that the simplest supersymmetric S'-matrix as above coincides 
with the S'-matrix coming from another integrable model. Consider, in fact, the S"-matrix 
in the solitonic sector of the ordinary Sine-Gordon model, given in full generality byQ 



S(P) = R{(3) 



sinh - 



sinh 



tv({3 — i-jv) 



sinh 



n({3 — i-jv) 







sinh 



sinh - 



sinh 








7T (p — in) 

e 








-1 





where 



R((3) = exp 



o 



°° dt sinh I (! - J 
t cosh | sinh p- 



sin 








(3.1 



71 



(3.9) 



2 The two-particle states are ordered as | SS), \ SS), \ SS), \ SS), where S denotes the soliton whereas 
S the anti-soliton. 
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The parameter £ is related to the coupling constant g of the Lagrangian of the model 



2 



1 7YI 

C = -(d^) 2 + —cos(gip) , (3.10) 



by£ 



2 



87T-9 2 ■ 

By comparing the S-matrix (|3.1|) with the one of the Sine-Gordon model, it is now 
easy to see that the two coincide for the particular value £ = 2ir, i.e. g 2 = 167r/3. This is 
a repulsive point of the Sine-Gordon model, no additional bound states are present at this 
value and therefore the full S'-matrix of the SG model is reduced to the one in (|3.8| ). Since 
the ordinary Sine-Gordon model can be regarded as a massive integrable deformation of 
the gaussian conformal model with an action 

A*=\j d 2 x(d»v) 2 , (3.11) 

(of central charge C — 1), the deforming operator [e %9Lp + e~ l9lp ) which leads to the Sine- 
Gordon model has at the point £ = 2tt the conformal dimension A = 2/3. The relevance 
of these observations as well as the consequences of the identity between the ^-matrices 
will be topics of discussion later in this paper. 

4 The SShG and Superconformal Models 

The aim of this section is to illustrate several properties of the Super Sinh-Gordon model 
as well as some features of the superconformal models and their deformations. 

4.1 Lagrangian of the SShG model 

In the euclidean space, the Super Sinh-Gordon model can be defined in terms of its action 



j <ui 

+ i— cosh A$ 
2 A 



(4.1) 



A = I d 2 zd 2 6 

where the covariant derivatives are defined as 

D = d g -6d z ; . , 

_ _ ' 4.2 

D = dg-9d g , 

and the superfield $(2:, z, 6, 9) has an expansion as 

z, 9, 9) = <p(z, z) + 9tjj(z, z) + 9lj>(z, z) + 99F(z, z) . (4.3) 

The integration on the 9 variables as well as the elimination of the auxiliary field !F(z, z) 
by means of its algebraic equation of motion leads to the Lagrangian 

1 ui 2 

£ = -(d z (pdz(p + ipd z ijj + ipdzifj) + -— sinh 2 \<p + inwpip cosh \<p . (4.4) 

2 2X 2 
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Of all the different ways of looking at the SSHG model, one of the most convenient is 
to consider it as a deformation of the superconformal model described by the action 

A,= l -f (d zV d- zV + tfdrf + i/>dd>) . (4.5) 

This superconformal model has central charge C = 3/2. At this poit, it is useful to briefly 
remind some properties of the superconformal models and their deformations. 



4.2 Superconformal models and their deformation 

For a generic superconformal model, the supersymmetric charges can be represented by 
the differential operators 

Q = d e + 9d z ; 

Q = dg + 6d- z . 

The analytic part of the stress-energy tensor T(z) and the current G(z) which generates 
the supersymmetry combine themselves into the analytic superfield 

W (z, 6) = G(z) + 6T(z) , (4.7) 

which is called the super stress-energy tensor. For the anti-analytic sector we have corre- 
spondingly W(z) = G(z) + 6T(z). These fields are mapped one into the other by means 
of the super-charges 



T(z) = {G(z),Q} , d z G=[T(z),Q] ; 
T(z) = {G(z),Q} , d- z G=[T(z),-Q] , 



(4.8) 



and their Operator Product Expansion reads 



C T(w) dT(w) 

T(z)G(w) = § GH + SGH + ... (49) 
2 [z — wY z — w 



G{z)G{w) = - 2 ° +2^ + - 
6{z — w) A (z — w) 



with analogous relations for the anti-analytic fields. As it is well known [p5|, reducible 
unitary representations of the N — 1 superconformal symmetry occurs for the discrete 
values of the central charge 

C = - - / 2 , . (4.10) 
2 m(m + 2) v ; 

At these values, realizations of the N = 1 superconformal algebra are given in terms of a 
finite number of superfields in the Neveu-Scwartz sector and a finite number of ordinary 
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conformal primary fields in the Ramond sector. Their conformal dimensions are given by 

^= [(ro+ g X" + 1'" 4 + ^ [1 - ( - 1 ^- (4 ' u) 

where p — q even corresponds to the primary Neveu-Schwartz superfields J\f(™>(z, 9) and 
p — q odd to the primary Ramond fields R^(z). These fields enter the so-called super- 
conformal minimal models SM. m . 

The Witten index Tr(— 1) F of the superconformal models can be computed by initially 
defining them on a cylinder P5[ : the Hamiltonian on the cylinder is given by H — 
Q 2 = L — C/24, where as usual C/24 is the Casimir energy on the cylinder and L = 
Y~ i dzzT(z). Considering that for any conformal state | a) with A > C/24 there is the 
companion state Q \ a) of opposite fermionic parity, their contributions cancel each other 
in Tr(— 1) F and therefore only the ground states with A = C/24 (which are not necessarily 
paired) enter the final expression of the Witten index. For the minimal models, there is 
a non-zero Witten index only for m even. Therefore the lowest superconformal minimal 
model with a non-zero Witten index is the one with m = 4, which has a central charge 
C = 1 and corresponds to the class of universality of the critical Ashkin- Teller model. 
The superconformal theory with C = 3/2 made of free bosonic and fermionic fields also 
has a non-zero Witten index, because an unpaired Ramond field R{z) is explicitly given 
by the spin field a(z) of Majorana fermion ip(z) with conformal dimension A = 1/16, i.e. 
by the magnetization operator of the Ising model. 

The above observations become important in the understanding the off-critical dy- 
namics relative to the deformation of the action of the superconformal minimal models 
SA4 m by means of the relevant supersymmetric Neveu-Schwartz operator Afl^\z,9). In 
fact, as shown in [26|, the massless Renormalization Group flow generated by such an op- 
erator preserves the Witten index. Therefore the long distance behaviour of the deformed 
SM., m minimal model - controlled by the action 

A„ + 7/ d 2 zd 2 9N[f , (4.12) 

(for a positive value of the coupling constant 7, with the usual conformal normalization 
of the superfield) - is ruled by the fixed point of the minimal model SM. m _2- Therefore 
the action (|4.12|) describes the RG flow 

SM m -> SM m -2 , (4.13) 
with a corresponding jump in steps of two of the central charge]], i.e. 

AC = C(m) - C(m - 2) . 

3 This is in contrast of what happens in the ordinary conformal minimal models, where the deformation 
of the conformal action by means of the operator ^3 induces a massless flow between two next neighbored 
minimal models. 
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Therefore, a cascade of massless flows which start from C = 3/2 and progress by all 
A/x 3 deformations of SA4 m met along the way must necessarily pass through the model 
with C — 1 in the second to last step rather than the lowest modefl with C = 7/10 
(Figure 1). We will see that this is indeed the scenario which is described by a specific 
analytical continuation of the coupling constant of the Super-Sinh-Gordon model, the 
so-called Roaming Models. 

For the time being, let us further discuss the deformation of the C = 3/2 super- 
conformal theory which leads to the SShG model. At the conformal point, the explicit 
realization of the component of the super stress-energy tensor are given by 

T(^) = -i[(^) 2 -^] ; 
G(z) = iipd z ip , 

and they satisfy the conservation laws d z T(z) = d z G(z) = 0. Once this super conformal 
model is deformed according to the Lagrangian (|4.4j) , the new conservation laws are given 
by 



d z T(z,z) = d z Q(z,z) ; 
d z G(z,z) = d z x(z,z) , 



(4.15) 



where 



Q(z,z) = — — sinh 2 Xip + imtjjtjj cosh Ay? , (4-16) 
IX 
j] j 

x(z,z) = —ip sinhAy? . 
A 

For the anti-analytic part of the super stress-energy tensor we have 

d z f(z,z)=d,Q(z,z) ; 
d z G(z, z) = d z x(z, z) , 

where Q(z, z) is as before and the other fields are given by 



(4.17) 



G{z,z) = -iipd z <p ; (4.18) 

TYl 

x{z, z) = —ip sinh \(p . 
A 

The operators Q(z,z), x(z,z) and x{ z i z) belong to the trace of the supersymmetric 
stress-energy tensor and they are related each other by 

Q(z,z) = {x(z 1 z),Q} , d zX (z,z) = [Q(z,z),Q] ; 

e(z,z) = {x(z,z),Q} , d,xM) = [e(z,z),g] , 



4 As well known, the lowest model with C — 7/10 corresponds to the class of universality of the 
Tricritical Ising Model. 
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where the charges of supersymmetry are expressed by 

Q = J G(z, z)dz + x(z, z)dz ; ^ 
Q = J G(z, z)dz + x(z, z)dz . 

In addition to the conservation laws ( |4.15| ) and ( |4.17| ), the SShG model possesses higher 
integrals of motion which were explicitly determined in p3| . Therefore its scattering 
processes are purely elastic and factorizable. and its two-body S-matrix is discussed in 
the next section. 



5 The S'-matrix of the SSHG 



The S'-matrix of the SSHG model has been determined in [24|. It is given by 



S(P) = Y((3) 



1 



2i sin net 
sinh (3 



cosh ^ 







cosh 2 



2i sin net 
sinh j3 








sinh ^ 



V 



o 



sinh ■ 



(5.1) 



where 



and the function U(f3) is given by 
U(J3) = exp 



sinh ~ 



sinh | + i sin na 



U(P,a) 



(5.2) 



dt sinh at sinh(l — a)t fit 

2~t sm — 

o t cosh 2 cosh t ix 



(5.3) 



The angle a is a positive quantity, related to the coupling constant A of the model by 

1 A 2 



a 



(5.4) 



This equation implies that the SShG is a quantum field theory invariant under the strong- 
weak duality 

47T 

A - T . (5.5) 

It is easy to see that such S'-matrix fulfills all constraints of section 2. The prefactor of 
the S-matrix admits the following infinite product representation 

1 



r(^£)r(^£ + i)r(| + i£)r(i + i£ 



°° r (k + a + \ 

xii 

k=Q 

oo 

><n 

k=0 



*&)r(* + 



a + i£)r*(k + i + i£ 



r(k + a + l-i£)r(k + 

T(k + a-i^)r(k + l- 



a 



Z7T 



(5.6) 



a 



r (* + a + l + i|L)r(* + 2-a + i£)ra (*+§ + *£ 
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from which one can explicitly sees that this S-matrix does not have poles into the physical 
strip. It is now interesting to analyse several analytic continuations of the above S'-matrix 
in the parameter a. 



5.1 The analytic continuation a — ► —a and Schoutens's model 



Under this analytic continuation, the SShG model goes into the Super Sine-Gordon model 
and the S'-matrix ( |5.1| ) describes in this case the scattering of the lowest breather states 
of the latter model. The explicit expression is given by 



S(P) = Y{(3) 



1 + 



2i sin not 
sinh /3 

sin ira 



cosh ^ 



cosh 1 





-1 + 



2i sin7ro 



sinh f3 

















sinh 2. 






1 



sinh ^ 



(5.7) 



where 



and 



sinh | 



sinh | — i sin ira 



U{(3) 



U{(3) = U{(3) 



I sinh | — i sin ira \ / sinh (3 + i sin(27r«) \ 
\ sinh | + i sin ira J \ sinh (3 — % sin(27r« J 



(5.? 



(5.9) 



Notice that at the particular value a = 7r/3 we can consistently truncate the theory at 
the Super Sine-Gordon breather sector only|. At this value, the pole at (3 = 27ri/3 of the 
S-matrix can be regarded as due to the bosonic and fermionic one-particle states | b(/3)) 
and | f{(3)). These particles are therefore bound states of themselves, in the channels 



bb -> b - 
ff^b 

for the bosonic particle b, and in the channels 



bf 
fb 



f 
f 



bb ; 
>// 



for the fermionic particle /. There is of course a price to pay for this truncation: this 

means that the residues of the S-matrix at these poles will be purely imaginary. Such a 

5 The general class of models considered by Schoutens is obtained by taking a — tt/(2N + 1), 
N = 1,2, .. . and they correspond to the supersymmetric deformation of the non-unitary minimal su- 
perconformal models with central charge C — — 3iV(4iV + 3)/(2iV + 2). For simplicity only the first is 
examined, the detailed discussion of the others requires the application of the bootstrap equations to 
their S'-matrix. 
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model therefore would be the supersymmetric analogous of the Yang-Lee model for the 



ordinary Sine-Gordon model p2| . It was considered originally by Schoutens |20| and it 
has been identified with the off-critical supersymmetric deformation of the non-unitary 
superconformal minimal model with central charge C = —21/4. The residues of the 
5-matrix are given by (Figure 2) 



1 bb 



1 // 

1 bf 
1 fb 



i V 3 v^3 k ; 
: % V3 k ; 

i \/3k , 



(5.10) 



where 



cxp 



°° dt sinh | sinh 



2 2t 
3 



t cosh 2 | cosh t 



0.7941001... 



(5.11) 



We will come back to this model for the discussion of its Form Factors. 



5.2 The Roaming Models 



The S matrix ( |5.1| ) of the SShG model has zeros in the physical strip located at a\ = 
ma and a 2 = 27r(l — a). By varying the coupling constant A, they move along the 
imaginary axis and they finally meet at the point «7r/2, at the self-dual value of the 
coupling constant A 2 = An. If we further increase the value of the coupling constant, they 
simply swap positions. But there is a more interesting possibility: as first proposed by 
Al. Zamolodchikov for the analogous case of the ordinary Sinh-Gordon model [p7| 



once 



the two zeros meet at in, they can enter the physical strip by taking complex values of 
the coupling constant (Figure 3). In this way, the location of the two zeros are given by 

a± = — ± ia . (5.12) 

From the analytic ^-matrix theory, the existence of complex zeros in the physical strip 
implies the presence of resonances in the system. By analysing the finite-size behaviour 
of the theory by means of the Thermodynamic Bethe Ansatz ||28|| , the interesting result 
is that the net effect of these resonances consists in an infinite cascade of massless Renor- 
malization Group flows generated by the Neveu-Schwartz fields J\f™ 3 and passing through 
all minimal superconformal model SM. m with non-zero Witten index (see Figure 1). As 
discussed in the previous section, the ending point of this infinite-nested RG flow should 
describe the deformation of the superconformal model SM.^. Is this really the case? 
By taking the limit a — > ioo into the S- matrix ( |5.1|) , it is easy to see that it reduces to 
the simplest supersymmetric S'-matrix analysed in Section 3, which in turn coincides with 
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the one of the Sine-Gordon model at £ = 2n. Since this S-matrix describes a massive 
deformation of the C = 1 model, in order to confirm the above roaming trajectory sce- 
nario the only thing that remains to check is the comparison of the anomalous dimension 
of deforming field. In the Sine-Gordon model at £ = 2ir, the anomalous dimension of the 
deforming field was determined to be A = 2/3, which is indeed the conformal dimension 
of the top component of the superfield jVip in the model SM.^ ! In the light of this result, 
it is now clear why in the roaming limit the value which is actually selected is £ = 2tt 
among all possible values of the coupling constant of Sine-Gordon model. 

In the next section we will see that the identity between the S-matrix of the two 
models also implies an identity between the Form Factors of the two theories. 



6 Form Factors of the Trace Operators of the SShG 
Model 

For integrable quantum field theories, the knowledge of the S'-matrix is very often the 
starting point for a complete solution of quantum field dynamics in terms of an explicit 
construction of the correlation functions of all fields of the theory. This result can be 
obtained by computing first the matrix elements of the operators on the asymptotic 
states (the so-called Form Factors) 0, [2] and then inserting them into the the spectral 
representation of the correlators. For instance, in the case of the two-point correlation 
function of a generic operator 0(z,z) we have 



Q{z,z) = (0 | O(z,z)O(0,0) | 0) = / da 2 p{a 2 ) K (ay/zl) , (6.1) 

Jo 

where K (x) is the usual Bessel function. The spectral density p(a 2 ) is given in this case 
by 

2 = g I d^h ^ a _ £ m cogh £ ginh x 

n=o J 2n 2n i 

(0 | 0(0,0) I A 1 ((3 1 )... A n {(3 n ))\ 2 . 



The Form-Factor approach has proved to be extremely successful for theories with scalar 
S'-matrix, leading to an explicit solution of models of statistical mechanics interest such as 
the Ising model [jl], [| [L4| , the Yang-Lee model |T3| or quantum field theories defined by a 



lagrangian, like the Sinh-Gordon model |7], ||, [13| . On the contrary, for theories with a non- 
scalar S-matrix the functional equations satisfied by the Form Factors are generally quite 
difficult to tackle and a part from the Sine-Gordon model or theories which can be brought 
back to it ]l], |2|, |J , there is presently no mathematical technique available for solving them 
in their full generality. Also in this case, however, the situation is not as impractible as 
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it might seem at first sight. The reason consists in the fast convergent behaviour of the 
spectral representation series which approximates the correlation functions with a high 



level of accuracy even if truncated at the first available matrix elements JTTj, [LJ, [U|, 16| [17 



18], In the light of this fact, in this section we will compute the lowest matrix elements 
of two of the most important operators of the theory, namely the trace operators Q(z, z), 
x(z, z) and x( z -> °f the supersymmetric stress-energy tensor of the SShG model. For 
the operator 9(0, 0), they are given byQ 

F*(fl = (0\e{0,0)\b(p 1 )Kh)) ; 



Ff f {p) = (0 | 9(0) | /(A)/(A)> 



(6.3) 



whereas for the operators x(0, 0) we have instead 



= (0 | x(0,0) | b{^)f{(5 2 )) ; 
Ff b (f3 1 ,(3 2 ) = (0\ X (0,0)\f(f3 1 )b((3 2 )) , 

(with an analogous result for the lowest Form Factors of the operator x(0,0)). Since 
the operators 9, \ and % are related each other by supersymmetry, as consequence of 
eqs. Q4.19| ) we have 



F b e b (P)= u;(e^Ff b + e^F b ) j 

Ff f (f3) = —To (e^ 2 F% - e^ 2 F§) ; (6.5) 
Ff Jfl) = - u ( e ~^ 2 F% - e~^ 2 F* 



It is therefore sufficient to compute the two-particle Form Factors of the operator Q(z, z) 
for determining those of x{ z i z) and x(z,z). Let us discuss the functional equations 
satisfied by Fg(J3) and Ff f (J3). 

The first set of equations (called the unitarity equations) rules the monodromy prop- 
erties of the matrix elements as dictated by the S'-matrix amplitudes 

F$(0) = S b b m F$(-0) + SUW Ff f {-(3) ; (6.6) 
Ff f {P) = S%(fi) F*{-ff) + S f / f (P) Ff f (-/3f) , 

where S bb (f3) is the scattering amplitude of two bosons into two bosons and similarly for 
the other amplitudes. 

The second set of equations (called crossing equations) express the locality of the 
operator Q(z, z) 

F 6 f(/3 + 2vrz)=F 6 f(-/3) 
if f (J3 + 2m)=Ff f { 



F? f (P + 2m)=Ff f (-(3) . (6 ' 7) 



5 They depend on the difference of rapidities (3 = fi\ — 02 since 6(0,0) is a scalar operator. 
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One may be inclined to solve the coupled monodromy-crossing equations by initially 
diagonalising the S-matrix However, this method does not work in this case for a 
series of reasons, both of mathematical and physical origin. To simplify the notation, let 
us use in the following 

a = sin rca ; 

s = sinh — ; 

2 ' 

c = cosh ^ . 

Observe that the eigenvalues of the S'-matrix of the SShG model in the F = 1 sector are 
given by 

C± = -^ ± V 1+ (c) ■ (6 ' 8) 
with the corresponding eigenvectors 

I C+> = 

I C-> = N. 

and the normalization constants given by 




From a mathematical point of view, the branch cuts present in the eigenvalues ( |6.8| ) make 
it impossible to find their exponential integral representations - a step which is usually 
rather crucial in obtaining the corresponding Form Factor [|IJ. A more serious aspect, 
however, is the fact that the eigenvectors Q6.9Q do not have any satisfactory properties 
under the crossing transformation (3 — > (3 + 2iri. From a physical point of view, the 
origin of all these troubles is the different scattering property of the channel bb — ► bb 
with respect to the channel // — > ff, which does not permit, in this case, to assign a 
reasonable physical meaning to the states which diagonalise the S'-matrix. 

The determination of the Form Factors of the operator Q(z,z) must pass through 
a different route. The way that we will proceed is to introduce two auxiliary functions 
F±(f$) by considering the scattering theory in the F = — 1 sector and to use them as 
building blocks for constructing the matrix elements F^{(3) and FfA(3). 



a 
c 

a 

c 



KPi)KP 2 )) +U- ^/l+(£) 2 j | /(&)/(&)) 

WTO) + ( i + I /(A)/(A)> 



(6.9) 
(6.10) 
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6.1 Auxiliary problem: two-particle FF in the F = — 1 sector 

Let us look for the eigenvalues of the ^-matrix in the F — — 1 sector. They are given by 

X + = S -^U((3) ; (6.11) 
s + ta 

A- = -U{p) , 
with the relative eigenvectors given by 

i a + > = -L (i b(/h)fm) + fwm))) ; (6-i2) 

I A^i(|/(AW 2 ))-W(ft)>) • 

The states in this sector are necessarily coupled to operators with a non-zero fermionic 
quantum number. For the purpose of obtaining the auxiliary functions F±(f3) which will 
be used as building blocks to construct F^(/3) and FfJfl), it is sufficient to consider the 
coupling of the eigenstates ( |6.12|) to a fictitious scalar operator A(0) but with a non-zero 



fermionic quantum number. Let us denote the corresponding matrix elements as 

F + (/3) = (0|A(0)|A + (/3 1 ,/3 2 )) ; 
F_(/3) = (0| A(0) | X-tfufo)) . 

Under the condition of unitarity, these matrix elements satisy the equations 

F + ((3)= S -^U(/3)F + (-(3) ; (6.14) 
s + la 

F-{0) = -U{p)F-(-0) . 

The crossing properties of the above matrix elements is more subtle. Let us consider, in 
fact, what happens to F + (f3) and to F_(f3) in the analytic continuation (3 — > (3 + 27ri, 
i.e. when the first particle in each of the two states entering the eigenvectors (|6.12 ) goes 



around the operator A(0). Since this operator has a fermionic quantum number, if the 
particle which goes around the operator is also a fermion we get an extra phase (— 1), 
otherwise nothing, and therefore 

(0 | A(0) (| m)fW)+ I /(&)&(&)» — "(0 I A(0) (| b((3 2 )f(f3 1 ))- | f(f3 2 )b(M)) 
(0 | A(0) (| m)W))- I &(&)/(&)» "(0 I A (0) (I b((3 2 )f(f3 1 ))+ | rwm))) 

(6.15) 

i.e. under crossing the two matrix elements mix each other, 

F + (/3 + 27ri) = -F_(-/3) ; 
F.((3 + 2 m ) = -F + (-f3) . 
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In order to solve the coupled system of crossing-unitarity equations ( |6.14| ) and ( |6.16| ), it 
is convenient to separate the problem into two steps. The first step consists in finding a 
function G(f3) which solves the monodromy equation involving only the function U(/3), 
with the usual crossing property, i.e. 

G((3) = U((3)G(-P) ; 
G((3 + 2m) =G(-P) 

The explicit expression of the function G(/3) can be found in Appendix A. 

The second step consists instead in finding two functions f+(P) and f-(0) which solve 
the functional equations 

s — ia 



(6.17) 



/+(£) 



s + ia 

/_(/?) = -f-(-P) ; (6-18) 

f+(P + 2iri) = -f-(-0) ; 
f_((3 + 2m) = -f + (-(3) . 

To this aim, let us write initially the integral representation of the eigenvalues A+(/3) 
s — ia 



s + ia 



cxp 



-Ai 



dt sinh at smhfl — a)t 3t 
sin — 

t COSh t 7T 



The minimal solutions of fl6.18|) are then given by 

P 



f+(j3) = cosh h -H + ((3) ; 
/_(/?)=* sinh ^tf_(/3) . 



(6.19) 



(6.20) 



where 



H + (f3)=exp 



H-{0) = exp 



00 dt sinh at sinh(l — a)t 

t cosh t sinh 2t 
00 (it sinh at sinh(l — a)t 



sm 



' (3-2m 



2n 



(6.21) 



• 2 /? 
sm — t 

2vr 



t cosh t sinh 2t 
The function H_(f3) admits the equivalent representation 

r (* + l + *£) r (* + 1 -f + *£)r(fc + i + f + 



^-09) = n 

fc=0 



Itt 



r(* + i + <£)r(* 



2 ^ Mtt 



r (fc + i + 



4-7T 



(6.22) 



or the mixed one, which is more useful for numerical calculation 



x exp 



JV-l 

n 

k=0 



1 + (7^+2) ) 







7r(4fc+4+2a) 



/3 



(4fc+6+2a) 



1 + 



iL_V 



(4^+6) y 

(it sinh at sinh(l — a)t 



jL 



o t 



cosh t sinh 2t 



(4fc+2+2a) 

(N + 1 - iVe -4 ' 



JL 



(4fc+4-2a) 



fe+1 



X 



P 



-wt in 2 _ t 
2ix 



(6.23) 
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Since the functions H±(0) are related each other by the equation H + ((3 + 2m) = H_((3), 
the corresponding infinite product or mixed representations of H + {j3) follow easily. For 
large values of /3, H + ((3) and H-(/3) go the same constant, given by 

lim H±{(3) = exp 

So, summarising the results of this section, the final expressions of the functions F + (/3) 
and F_(/3) are given by 

F+(P) = G{p) /+(/?) ; 
F-(J3) = G(J3) /_(/?) , 

with given in eq.( |A.2|) of Appendix A and f±(P) given in eq. ( |6.20| ). 



00 cftsinhat sinh(l — a)t 
t cosh £ sinh It 



(6.24) 



6.2 Two-particle FF of the operator G 



Let us look for the two-particle Form Factors of the trace of the stress-energy tensor 
as linear combination of the two functions F + (f3) and F_{(5) above determined, i.e. 



F bb (P) = A(P)F + (P) + B(P)F.(P) ; 
F ff ((3) = C((3)F + ((3)+V((3)F.(f3) . 



(6.26) 



Let us now plug them into eq.( |6.6| ). By using the monodromy equations satisfied by the 
F±{(3) and by comparing the terms in front of each of these functions, we obtain the 
following equations 



A{(3) 

By taking the free limit a - 



s — la 

s 

s + ia 
s 

s — ia 
s 

s + ia 



1 - 



ia 

sc 
ia 



A{-(3)--C{-(3) 



sc 
ia 



l--)B{-0)--V{-0) 



(6.27) 



\ sc) c 



1 a 



sc 



in the above equations, it is easy to see that 

A(P) = A(-(3) ; 
B{(3) = -B(-P) ; 
C(f3) = -C{-(5) ; 
V(f3) = V{-(5) . 



(6.28) 



Let us assume that the above parity properties of the functions A, . . . , T> are still valid for 
a non-zero value of the coupling constant a (after all, the coupling constant dependence 
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of the matrix elements F bb and Fff should be already included into the functions F± (ft)). 
Under this hypothesis, eqs. (|6.27|) provide the relationships 



C{p) = -i — A{P) ; 



V(ft) 



s 

c + 1 



(6.29) 



B{ft) 



Hence, at this stage we have for the two-particle Form Factors of G(x) 



F bb ({3) = AF + + BF 
c- 1 



(6.30) 



Fff(ft) 



c+1 
AF + -—BF^ 



In order to determine the two remaining functions A(ft) and B(/3), let us consider once 
again the case a — > and let us impose the condition that in this limit the two-particle 
Form Factors reduce to their free limit 



F bb {(3) : 

Fff(ft) 



lirvrir ; 
— 27ri sinh | 



(6.31) 



From this matching, the functions A and B are uniquely determined to be 

A(ft) = 2ir cosh ^ , 
B{f3) = 2m sinh ^ . 



(6.32) 



For a generic value of the coupling constant, the correct expressions of the two-particle 
Form Factors of the operator O(0) can be otained by imposing their normalization F bb (i7r) = 
Fff{m) = 27rm 2 and their final form are given by 



where 



Mft) 



F b ^P) = 2^S^\ 
F bb (m) 

Ff f {3) = 2nm 2 



cosh 2 ^ H+(J3) - sinh 2 f H-(/3) 



(6.33) 



4 



G{ft) 



(6.34) 



ft 



Fff(ft) = sinh | [H + (P) + H4P)] G(ft) 



Notice that for large values of ft, F b ®((3) tends to a constant whereas Ff^{(3) ~ e 
both behaviour in agreement with Weinberg's power counting theorem of the Feynman 
diagrams. 
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While we postpone non trivial checks of the validity of ( p. 33 ) to the next sections, let 
us use the Form Factors ( |6.33j ) to estimate the correlation function C(r) = (0(r)0(O)) 
by means of formulas (|6.1|) and fl6.2| ). In the free limit, the correlator is simply expressed 
in terms of Bessel functions, 

C{r)=m i (Kl(mr) + Kl(mr)) . (6.35) 

For a finite value of a, a numerical integration of ( |6.1| ) produces the graphs shown in 
Figure 4. As it was expected, in the ultraviolet limit the curve relative to a finite value 
of a is stepest than the curve relative to the free case whereas it decreases slower at large 
values of mr. This curve is expected to correctly capture the long distance behaviour 
of the correlator and to provide a reasonable estimate of their short distance singularity. 
However, for the exact estimation of the power law singularity at the origin one would of 
course need the knowledge of all higher particle Form Factors. 



7 C-theorem Sum Rule 

As mentioned in Section 4.1, the SShG model can be seen as a massive deformation of 
the superconformal model with the central charge C — |. While this fixed point rules 
the ultraviolet properties of the model, its large distance behaviour is controlled by a 
purely massive theory with C — 0. The variation of the central charge in this RG flow is 
dictated by the C-theorem of Zamolodchikov [30(1, which we will discuss in more details 



in the next section in relation with its fermionic formulation. In this section, we are 
concerning with the integral version of the C-theorem in order to have non trivial 



checks of the validity of the Form Factors (|6.33|) . In this formulation of the c-theorem, 
the variation of the central charge AC satisfies the sum rule 

AC=—jd 2 x\x\ 2 (0\ 0(x)0(O) | 0) conn = j dfic(fi) , (7.1) 

where c(/z) is given by 



6 1 

7T 2 jJ, 

,2 



c(/i) = --ImG(p 2 = V) 



G(p 2 ) = / d 2 xe-* px (0 | 0(x)0(O) | 0} conn . (7.2) 



Inserting a complete set of in-state into ([7.2|), the spectral function c(/i) can be expressed 
as a sum on the FF's 

c(//) = ^ ^2 I ' " " 17 — ^(A* - ^mcoshPi) 5(^msinh/3j) x (7.3) 

^ n =l J 271 Z?l i i 



(0 | 0(0,0) | A x {Px)...AJfi n )) 



|2 
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Since the term | x | 2 present in ( |7. 1| ) suppresses the ultraviolet singularity of the two-point 
correlator of 6, the sum rule Q7.1J) is expected to be saturated by the first terms of the 
series f [7.3|) . For the SShG model the first approximation to the sum rule ( |7. 1|) is given by 
the contributions of the two-particle states 

The numerical data relative to the above integral for different values of the coupling 
constant a is reported in Table 1. They are remarkably close to the theoretical value 
AC = |, even for the largest possible value of the coupling constant, which is the self- 
dual point A = \fAn. In addition to this satisfactory check, a more interesting result 
is obtained by analysing the application of the c-theorem sum rule to models which are 
obtained as analytic continuation of the SShG. 

7.1 C-theorem Sum Rule for the Roaming Model 

As we have seen in Section 5.2, by taking the analytic continuation 

a -> ^ + ia (7.5) 

and the limit «o ~~ > °°> the 5-matrix of the SShG model is mapped into the 5-matrix of 
the ordinary Sine-Gordon model at £ = 2n. This identity between the two S'-matrices 
is expected to extend to other properties of the two theories as well. There are however 
some subtleties in this mapping. 

One subtlety is that the right hand side of the sum rule (|7.1| ) is a pure number (AC = | 
for the SShG model) - a quantity which therefore seems to be completly insensitive to the 
variation of the coupling constant of the model. On the other hand, if the SShG model 
is mapped onto the Sine-Gordon model in the limit a ~~ > 00 of the analytic continuation 
([7.5p, the sum rule should instead jump discontinously from the value AC = | to AC = 1, 
simply because the latter model is a massive deformation of the conformal theory with 
C — 1. This is indeed what happens^], the reason being the presence of an additional 
energy scale brought by the parameter a an d the non-uniform convergence of the spectral 
series (|7.3|). 

In the case of the SShG model, there is an additional point that however deserves 
to be carefully checked. Note that in the SShG model, the lowest approximation of the 
sum rule, eq. ( |7.4|) , is expressed as sum of the moduli square of both F b f and Ff*. In the 

7 We will not discuss further this statement since a detailed discussion of this aspect can be found in 
j§] , where the analogous situation of the roaming limit mapping between the Sinh-Gordon model and the 
Ising model is analysed. 
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roaming limit, the original states | bb) and | //) of the Sinh-Gordon model are mapped 
respectively onto the states | SS) and | SS) of the Sine-Gordon model. On the other 
hand, the trace G of the stress-energy tensor of the Sine-Gordon model only couples to 
the symmetric combination of the soliton states 



+) = -j= (| SS) + | SS) 



(7.6) 



and not to each of them individually! How this apparent discrepancy can be settled? 

First let us consider how the two-particle Form Factor of the operator G is computed 
within the Sine-Gordon theory for a generic value of £ [f| . The symmetry responsable for 
the coupling of the operator G to the combination (|7.6|) is the invariance of Q under the 
charge conjugation. Hence, to compute the two-particle Form Factor 



F° G ((3) = (0 | 0(0) | +) , 



(7.7) 



one needs to employ the corresponding eigenvalues of the S'-matrix 

c _ sinh §((3 + in) A 

and solve the unitarity and crossing equations 



given by 



(7.8) 



Ff G (f3 + 2ni) = Ff G (-/3) . 



(7.9) 



The solution of these equation gives us the two-particle Form Factor of of the Sine- 
Gordon model 

sinh (3 



with 



G(/U)=exp 



V2n 2 m 2 



sinh — in) 



dt sinh f (l - |) 
1 1 t cosh ~ sinh |^ sinh t 



sin 



in — 
2tt 



(7.10) 



(7.11) 



Observe that for £ = 2n, the above expression (|7.10|) can be written equivalently as 



Ff G (/3) 



V2nm 2 . , /3 
— smh — 

G(i7r,27r) 2 



sinh — + i cosh — 
4 4 



G(/?,2tt) . 



(7.12) 



Once inserted into the sum rule (ffll) , the result is AC^ 2 ) = 0.9924..., i.e. a saturation 
within few percent of the exact value AC = 1 relative to this case. 

Let us consider now the corresponding FF of the SShG model in the roaming limit, 



eq. (|6.25| ). For cto — > oo, the function G(f3) reduces precisely to G((3,2n) whereas for 
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H±(fl) we have 



H + (f3)=ex V 
H_({3) = exp 



o t 



sinh 2t 

2 



(7.13) 



00 dt sin 2 £-£ 
o t sinh 2t 



By using the formula 



00 dt 



at 



sin 



In 



1 + 



4x 2 



/o t 2 
and the infinite product representation 

oo / 

coshx = TT 1 + . , , 
it is easy to see that they can be simply expressed as 



a 



H+(J3) 



-i sinh 



H-(J3) = cosh 







(7.14) 



4 ' v / 4 

Hence, in the roaming limit the two matrix elements F^{(3) and F^{(3) become equal, 
with their common value given by 



V2nm 2 , (3 
■ — sinh — 

G(i7r,27r) 2 



sinh ^- + i cosh ^ 
4 4 



(7.15) 



which coincides with the one of the Sine-Gordon model, eq. ( 7.12j) . Therefore their con- 
tribution to the c-theorem sum rule is precisely the same as the Sine-Gordon model (in 
the case of the SShG model, eq. (|7.4|), there is in fact a factor | with respect to the 
Sine-Gordon model, which is cancelled by the equality of the two Form Factors _F 6 ® and 



7.2 C-theorem Sum Rule for Schoutens's Model 

In the analytic continuation a — > —a, the S'-matrix develops a pole in the physical strip 
located at (3 — luai. As discussed in Section 5.1, the scattering theory for a — ^ admits a 
consistent interpretation in terms of a multiplet made of a boson and a fermion, which are 
bound states of themselves. For this model, the operator 0(z, z) has also a one-particle 
Form Factor F® in virtue of the graph of Figure 5. 

In order to discuss the analytic continuation of the Form Factors F^{(3) and FfAfl) 
for a —> —a and to calculate the matrix element F®, it is convenient to consider initially 
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a different representation of the functions f±{(3). This is given by 



/-(/*) 



sinh^Z+03) ; 
-sinh^Z_(/3) 



(7.16) 



with 



Z+(/3) = exp 

Z-W) 



exp 



dtcosh(l - 2a)t . 2 f{3 -2m 
t cosh t sinh 2t " \ 2n 



dt coshfl — 2a)t 3 

— f - sin 2 — i 

i cosh t sinh 2£ 27r 



(7.17) 



Their equivalence to the previous ones, eq. (|6.20|) is easily established by comparing their 
infinite product representations which for the function Z_(/3) is given by 

2 



Z_{3) 



n 

k=0 



T k 



2 ^ 4 4vr 



r fife 



1 ^ 2 ^4tt 



(7.18) 



while the one of Z + (3) is simply given by the identity Z + (8 + 2ni) = Z_(3). By us- 
ing eq. ([7.1 8|) and the analogous for Z + (8), it is simply to see that under the analytic 
continuation a — > —a, the functions /±(/3) change as follows 

2tt 2 



/-(/?) 



/+(/?) 
/-(/?) 



cosh £ — cos na 

ZZ\3) • 



(7.19) 



2 

27T 2 



cosh | + cos -na 



The analytic continuation of the two functions present therefore distinct physical prop- 
erties: in fact, while the function /+(/?) develops a pole relative to the bound state, the 
other /_(/?) does not have any singularity in the physical strip. 

To compute the two-particle FF of G is not sufficient, however, to substitute the above 
functions f±(3) into eq. ( |6.33| ). In fact, we must also take into account the extra terms 
present in U(3), eq. (|5.9| ). This is a simple task, though, because we only have to introduce 
the function K(8) = Kx(8)K 2 (8), where K\(8) and K 2 {8) are solutions of the unitarity 
and crossing equations 



K x {8) 



I sinh | — % sin na " 
\ sinh 77 + i sin ira / 



/ sinh 3 + i sin(27ro;) \ 
y sinh 8 — i sin(27ro; J 



(7.20) 



K 2 (8) 



K 2 {-8) 



K 1 (P + 2m) = K 2 (-P) . 
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Their mixed representation expressions are given by 



N-l 

km = n 

fc=0 



1 + 



in — ft 



1 + 



in — (3 



X 



x exp 



/ 

Jo 



dt cosh(l — 2a)t 



_ 4 Zl v- "V- e ~4nt gin 2 ll±t 



(Z7T-/3). 



and 



JV-l 

^ 2 (/3) = n 

fc=0 



i sinh 2t 



in— (5 
2vr(fc+i + f ) 



27T 



(7.21) 



17T-/3 



2vr(fc+l-§) 



x exp 



27r(fc+ 

00 cosh |(1 - 2a) 
o t sinh t cosh | 



iffy) ) (i + (i^F 



fc+i 



X 



(TV + 1 — iVe- 2 ')e- 2n * sin 2 ^ ^ 

27T 



(7.22) 



So, finally the two-particle FF of the operator Q(z,z) of the Schoutens's model are 
given by 



n&("0 

F° (/?) = 27rm 2 



(7.23) 



F }f(. ilx ) 



where 



F bb {(3) = sinh 







, P Z +\P) - , 

cosh — -; - 1 r- + % sinh 

4 (cosh | — COS7TQ:] 



F /7 (/3) = sinh 



P 



COS7TOJ 



P Z-\P) 

4 (cosh | + cos7ra 



sinh j -. + — r- - % cosh j - ( ^ — — 

4 (cosh I — cos7ra] 4 (cosh | + 



cos na 



G{P)K{(3) ; 
G(/3)X(/3) . 



The one-particle Form Factor F b can now be equivalently obtained either from the residue 
equation on F b %((3) or on Fff 



-i lim f/3 



3 ) 

2lxi \ ire 



(7.24) 



Ffm=T) f F® , 



with the result 



-1.6719(3) i 



(7.25) 
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The series of the sum rule has alternating sign, with the first contribution given by the 
one-particle Form Factor 

AC {1) = — (F®) 2 = -5.3387(4) (7.26) 

7T 



21 
' 4 

including the positive contribution of the two-particle FF, computed numerically 



This quantity differs for a 1.6% from the theoretical value AC = — = —5.25. By also 



AC (2) = 0.09050(8) , (7.27) 

the estimate of the central charge of the model further improves, C = —5.2482(4), with a 
difference from the exact value of just 0.033%. 

The Form Factors above determined for the trace operator Q(x) can be used to esti- 
mate its two-point correlation function. The graph of this function is shown in Figure 6: 
note that this function diverges at the origin, in agreement with the positive conformal 
dimension A = 1/4 of this operator, but it presents a non-monotonous behaviour for the 
alternating sign of its spectral series. 

8 Fermionic Formulation of the C-theorem 



In the OPE Q4.9|) of the fields T(z) and G{z) which are responsible for a superconfor- 
mal symmetry, the central charge C appears both in the short-distance singularity of 
T(z)T(w) and G(z)G(w). Therefore, going away from criticality by means of operators 
which preserve the supersymmetry of the critical point (but non necessarely its integra- 
bility), it should be possible to formulate a C-theorem for unitarity theories by looking 
at the fermionic sector. This is indeed the proved in this section. 

Let us consider the correlation functions of the operators G(z, z) and xi z i %) away 
from criticality. By taking into account the scaling dimensions of the operators, these 
correlators can be parameterised as 

H(zz) 

(G(z,z)G(0,0)) = 



zr 
Kzz) 

(G(z,z) X (0,0)) = ; (8.1) 

z z z 

( X (z,z)x(0,0)) = , 

ZZ A 

where the functions H(zz), I(zz) and L(zz) depend on the scalar variable x = zz. In 
virtue of the conservation law 

d g G(z,z)=d zX {z,z) , (8.2) 
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the above functions satisfy the following differential equations 

H - I = -21 ; 

. . 8.3 

I — L = I — L , v ' 

where H = zzj^ and similarly for the other functions. For the combination F = H+I — L 
we have therefore 

F = -2L . (8.4) 
At the critical points, the quantity F is proportional to the central charge C, since 

H = I C > (8.5) 
I=L=0 . V 1 

In order to conclude that the function F decreases along the Renormalization Group flow 
and to derive the corresponding sum-rule, it is necessary to analyse the positivity of the 
function L(zz), which is given by 

L(zz) = (zz)z(x(z,z) X (0,0)) . (8.6) 

Although it seems a-priori difficult to prove the positivity of the above function for the 
explicit presence of the complex term z, it will become evident that this term is precisely 
required for the proof. 

Let us start our analysis from the constraints provided by supersymmetry. By using 
eq. (|4.19| ), we have 

(Q(z, z)8(0, 0)> = ({ X (z, z), Q}{ X (0, 0), Q}) = (8.7) 
= (x(z,z)Q 2 X (0,0)) = (x(z,z)P +X (0,0)) • 

The above equations implies that the spectral function of the correlator (Q(z, z)Q(0, 0)) 
can be obtained from the one of (x(z, ^)x(0, 0)) by multiplying each term of the multi- 
particle expansion for the right momentum of the intermediate state. This is a positive 
quantity, since for a cluster of n-particles is given by 

P + = m (e pl + e^ 2 + ■ ■ ■ e Pn ) . (8.8) 

The positivity of the spectral series of (Q(z, z)Q(0, 0) for unitarity theories implies there- 
fore the positivity of the spectral series of (x{z, ^)%(0, 0)) as well. 

Let us consider now the behaviour of the correlator (x(z,z) x{0,0)) as a function of 
the variables z and z. Since (|8.7|) can be written equivalently as 



(Q(z, z)G(0, 0)) = ( X (z, z) |£ x (0, 0)) , (8.9) 
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where M 2 is the invariant mass square of the cluster of the intermediate particles, this 
implies that (x(z, z)x(0,0)) can be expressed as 

(X(MM0,0)) = f^P-X(p 2 )e-^ p+z+p -~ z) = (8.10) 

J ATI 

The spectral density x(p 2 ) is given in terms of the matrix elements of the operator x once 
we have factorised the left momentum of the intermediate states. According to the above 
discussion this is a positive quantity which depends only on the invariant mass square of 
the cluster of the intermediate particles. By using the identity 

POO 

2\ _ / J„2 r/„ 2 „2w„2\ 



X (p 2 )= / da 2 S(p 2 - a 2 ) X (a 2 ) , 
Jo 

eq. ( |8.10| ) can be written then as 

[ f) roo 

(x(z,z)x(0,0)) = —t- da 2 X (a 2 ) Ko(aV^) . (8.11) 



7T 8z Jo 

Since 

dK 

^ = ~ Kl{x) ' 
by taking the derivative of the above expression we have 



1 fz 



(x(z,z) X (0,0)) = ^J- I da 2 ax{a 2 )K x {a^JTz) . (8.12) 
2tt V z Jo 



and therefore for the function L we have 



L(zz) = ^ — / da 2 ax(a 2 ) K^ay/zz) . (8.13) 



2tt Jo 

Eq. ( |8.13| ) explicitly manifests the positivity of this function, Q.E.D. 

The relative sum rule is easily obtained. In fact, since at the fixed points H = |C, we 



have 



Since 



C m -C fm = ^Jd(zz)z( x (z,z)x(0 1 0)}= (8.14) 

= — / da 2 a X (a 2 ) / dRR 2 K 1 (aR) . 
Itx Jo Jo 

r 2 
/ dRR 2 K l (aR) = — , 
J a 



eq. (|8.14j ) can be written as 



3 r°° da 2 , 9 . . 
AC = - / — X (a 2 ) . (8.15) 



7T JO GT 
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Let us consider now the two-particle contributions of the spectral function x(a 2 ). Since 

= ^e~^F®{{3)+uje-^Ff f {{3)) ; (8.16) 



we have 



F e 12 , I i2 



4 cosh 2 | 



(8.17) 



To obtain the relative contribution for x(a 2 ) we have just to disregard the factor 

( e -A/2 + e -A/2) 

relative to the left component of the momentum, so that 



2n 



1 



5 (a — m cosh (3% — m cosh fy) 5{m sinh (3i + m sinh j3%) x 



4 cosh 2 



(8.18) 



-5(a - 2m cosh (3) \ F%(2(3) | 2 + | Ff f (2/3) 



167r J-oo cosh 3 /? 
Inserting this expression into eq. ( |8.15| ), we finally have 



AC® 



d(3 



7r 2 m 4 Jo cosh (3 



Ff f (2P) 



ff 



(8.19) 



i.e. the sum rule relative to the fermionic case employs the same integrand as the bosonic 
case and therefore provides the same results as before. 



9 Summary and Conclusions 

In this paper we have analysed several aspects of supersymmetric models from the point 
of view of the S-matrix approach. The simplest supersymmetric scattering theory dis- 
cussed in Section 2 has a natural interpretation and identification in terms of the first 
model described by the Roaming Series. We have determined the lowest Form Factors of 
the trace operators of the super stress-energy tensor and we have discussed the novelty 
present in this calculation. Several checks of their validity have been presented in Section 
7, in particular a remarkable saturation of the c-theorem sum rule obtained for the sim- 
plest Schoutens's model. We have finally obtained a fermionic version of the c-theorem 
by employing the supersymmetric properties of the correlators. It would be extremely 
interesting to explore further the properties of these supersymmetric theories and obtain 
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closed expressions for Form Factors of other operators and with an arbitrary number of 
external particles. 
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Appendix A 



In this appendix we collect some useful formulas of the function G(j3) used in Section 6. 
The function G(/3) is the minimal solution of the equations 



G(i3) = U(l3)G(-l3) ; 
G(in -(3) = G(in + (5) , 

where U(f3) is the function entering the S'-matrix of the SShG model, eq. 
expression is given by 



G(l3) = exp 



dt sinh at sinhfl — a)t 9 (in — 3)t 

■ sin 2 



o t cosh 2 | sinh t cosh t 



271 



For large value of /3, this function tends to a constant given by 

_,_ N 1 f°° dt sinhat sinhfl — a)t 

hm G{p) = exp — - / 27 

P-*°c 2 Jo t cosh 2 sinh t cosh t 



(A.l) 
Its explicit 

(A.2) 
(A.3) 



For the purposes of numerical calculation, it can be written in a more convenient way as 
follows. First of all, notice that the function U(/3) can be factorised as 



where 



U 2 (p) = exp 



—i 



—i 



U{(3)=U 1 {(3)U 2 {(5) 



dt (1 - cosh 0(1 + coshfl - 2a)t dt 

sin — 



o t 



sinh t 



71 



(A.4) 



(A.5) 



dt cosh(l - 2a)t (3t 
sin — 

t COSht 71 



Correspondingly, G can be factorised as 

G(f3)=G 1 ((3)G 2 (f3) 



where 



N-l 

gm = n 

fc=0 



pl((3,0)pl((3A-2a) 
f k (J3,l)p(J3,2a)p k (/3,2-2a) 



(k + l)(k + 2) 
4 



(A.6) 



(A.7) 



x exp 



dt (1 - cosht)(l + cosh(l - 2a)t)d N (t) ^ 2 (in - l3)t 



o t 



sinh 6 1 
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and 



JV-l 

gm = n 

k=0 



{in -13) 



(Z7T-/3) 



x exp 



\(4k + 1 + 2a)n / 
dt cosh(l - 2a)te~ mt 



(4k + 3-2a)n j 



x (A.i 



o t 



cosh t sinh t 



sin 



(in - (3)t 
2n 
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with 



and 



(wr - g) x 
(2A; + 3 + x)7r 



(«r - g) 
(2fc + 3 - x) 



7T 



1 r 



d^(t) = | [(JV + 1)(JV + 2) - 2iV(iV + 2)e" 2 * + iV(iV + l)e" 4 * 



-2iVt 



Appendix B 

In this appendix we briefly describe the reason of the inconstency of the Form Factors 
found by Ahn [^J] for F 6 f and F^- of the SShG model, referring the reader to his original 
paper for all notations used by this author^ 

i.e. 



The Form Factors found by Ahn are those of formula (4.24) of 24 



-F_ mm ) , /3 

cosh — 

2 2 



Ff f {f3) = 2m' 



(jpmin _|_ pmin\ a 

' n 1 Tr-t Vt 



where 



= exp 



F_(p) = exp 



sinh — 
2 



o t sinht 



2tt 



dtf-(t) . 2 (iTT-^t 

sin 



o t sinht 



2tt 



and 



f±(t) 



[1 - cosht)(l + cosh((l - 2 | a \)t)) . cosh(l - 2 | a \)t) 



± 



(B.l) 



(B.2) 
(B.3) 



(B.4) 



shuV £ cosh t 

By construction, the functions F±(/3) considered by Ahn satisfied the unitarity equations 
relative to the eigenvalues of the S"-matrix in the F = — 1 sector and the crossing equations 
in the usual way, namely 



F + (P + 2m) = F + (-P) ; 
F_(/3 + 27ri) = F-(-/3) . 



(B.5) 



As discussed in Section 6, the crossing properties of the Form Factors in the F = — 1 on 
the other hand are quite different. This has an important consequence on the validity of 



8 They slightly differ from ours. The reader should also be aware that there are several misprints in 
the paper pi . 
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( |B.1| ). Consider in fact the transformation (3 — > (3 + 2ni in the first of (BTT). It is easy to 
see that the Form Factor of the trace operator 6 does not satisfy in this case the equation 



F 6 f(/3 + 27rz)=F 6 f(-/3) , (B.6) 

expressing the locality of this operator. 

The problem becomes evident as one tries to apply the c-theorem sum rule. If this 
is applied by using the above formulas and equation (|7.4j) , one gets the results reported 
in Table 2 which are definitely in disagreement with the theoretical result C = 3/2. The 
results reported in the original Table 1 of the paper by Ahn were actually obtained by 
means of an unjustified step at this stage of his calculations, namely by plugging negative 
values for the quantity | a | which appears in all formulas of Ahn's paper (see, for instance 
eq. ( |B.4|) above). In this case, although the first values reported by Ahn seem in reasonable 
agreement with the c-theorem sum rule, if that calculation had been pursued for higher 
values of the coupling constant, the results would have been quite unsatisfactory, as shown 
by the unreasonable small values obtained for the central charge for certain cases and a 
weird increasing behaviour around | a \— —0.5, see Table 3. 
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Table Caption 



Table 1 . The two-particle contribution to the sum rule of the c-theorem for the SShG 
model. 

Table 2 . The two-particle contribution to the c-theorem sum rule relative to the Ahn's 
Form Factors. 

Table 3 . The two-particle contribution to the c-theorem sum rule relative to the Ahn's 
Form Factors for negative values of | a \ . 



38 



a 


4tt 


Ac< 2 ' 


precision% 


1 

100 

3 
100 

1 

20 


1 

99 
3 

97 
1 

19 


1.49968 
1.49741 
1.49349 


0.0213 
0.1726 
0.4340 


1 

10 


1 

9 


1.47955 


1.3633 


3 

20 


3 
17 


1.46333 


2.4446 


1 
5 


1 

4 


1.44742 


3.5053 


3 

10 


3 
7 


1.42109 


5.2606 


2 
5 
1 
2 


2 
3 

1 


1.40480 
1.39935 


6.3466 
6.7100 



Table 1 
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ot 



1 


1 


100 


99 


3 


3 


100 


97 


1 


1 


20 


19 


1 


1 


10 


9 


3 


3 


20 


17 


1 


1 


5 


4 


3 


3 


10 


7 


2 


2 


5 


3 


1 
2 


1 



Ac< 2 ' 



1 


51481 


1 


54328 


1 


57019 


1 


63080 


1 


68227 


1 


72524 


1 


78833 


1 


82442 


1 


83616 



Table 2 
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a 


Ac( 2 ) 


-0.001 


1.4985 


-0.005 


1.49246 


-0.01 


1.48485 


-0.02 


1.46951 


-0.03 


1.45413 


-0.05 


1.42369 


-0.1 


1.35190 


-0.15 


1.28762 


-0.20 


1.23038 


-0.33 


1.14079 


-0.40. 


1.14552 


-0.50. 


1.58825 



Table 3 
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Figure Captions 

Figure 1 . Renormalization Group Flows described by the Roaming Models. 

Figure 2 . Bound states and residues in some amplitudes of the Schoutens's S-matrix. 

Figure 3 . Complex zeros of the Roaming Models. 

Figure 4 . Logarithm of the correlation function m _4 (O(:r)O(0)) versus mr for a = 
(full line) and for a = 0.5 (dashed line). 

Figure 5 . Bound state recursive equation for the form factor T 2 - 

Figure 6 . Correlation function m _4 (6(x)6(0)) versus mr of the first Schoutens's model, 
obtained with the first two terms of the spectral series. 
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Figure 4 
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Figure 6 
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